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We present an investigation into effects exhibited by the two-frequency kicked rotor. Experiments 
were performed and in addition quantum and classical dynamics were simulated and compared 
with the experimental results. The experiments involved pulsing the optical standing wave with 
two pulsing periods of differing frequencies and variable initial phase offset. The ratio of pulsing 
periods was sampled for rational and irrational values for different experimental runs. In this paper 
we present these results and examine the measured momentum distributions for the cause of any 
structures that are seen in the energy as the initial phase offset is changed. Irrational ratios exhibit 
no significant quantum effects, whereas rational ratios show dynamical localisation (DL) for certain 
values of the initial phase. However, most of the observed structure is found to be due to classical 
effects, in particular KAM boundaries, and is therefore not of uniquely quantum origin. 

PACS numbers: 



I. INTRODUCTION 

One of the signatures of quantum behaviour in the 
atom optics kicked rotor (AOKR) is dynamical localisa- 
tion (DL). For a classically chaotic phase space, the en- 
ergy of the classical kicked rotor increases approximately 
linearly with time and grows indefinitely. However, in the 
corresponding quantum system, the initial energy diffu- 
sion is approximately linear, but after some characteristic 
(localisation) time the energy growth is curtailed 0, • 
This halting of energy growth is due to quantum me- 
chanical interference and is expected to be very sensitive 
to changes in periodicity, as the quantum correlations 
are sensitively dependent on the time between successive 
kicks 0. 

It has been shown experimentally that DL is sensitive 
to decoherence in the form of spontaneous emission 0,H|. 
Previous work has also shown that in the presence of a 
small amount of period noise, DL is no longer observed 
000] Experimental results in 0, where two-frequency 
driving was employed, suggest that DL may only be ob- 
served when the ratio of the driving frequencies is ra- 
tional. These same results also exhibit sub-Fourier reso- 
nances [io|. 

In the experiment performed in [j| the ratio, r', of the 
pulsing frequencies of the two pulse-trains was changed 
and the number of zero velocity atoms measured. The 
time between the first pulse of each pulse-train, relative 
to the period of the second pulse-train (the phase), was 
set to 52° (where 360° corresponds to one period of the 
second pulse train) to avoid pulse overlaps. At highly 
rational values of r' (1/2, 1, 2, etc.) a large increase was 
observed in the number of zero velocity atoms. This was 
attributed to DL, as the same effect was not observed in 
simulations of the classical (S-kicked system. 

If pulse overlaps were able to be taken into account, the 
experiment would not be limited to any particular initial 
phase and the effective periodicity could be altered by 
changing that phase. In particular, investigating phase 



dependence would ensure that the results in [9j are not 
merely due to a particular choice of parameters, and that 
other effects do not come into play. 

In this paper we again consider the case where the 
driving may be periodic or quasi-periodic. We add a sec- 
ond pulsing frequency (period T2) to our primary pulsing 
frequency (period T\), as in Since DL is extremely 
sensitive to changes in periodicity, very small changes in 
the initial phase offset, ipo, should affect the amount of 
DL observable on the finite experimental timescale. This 
effect should occur for commensurate frequencies (as DL 
occurs for pulse-trains with sufficient periodicity), but 
not for incommensurate frequencies. 

The ratio r = T2/T1 defines the period of the second 
pulse-train. The rational or irrational nature of r should 
determine whether or not DL is observable on a particu- 
lar time scale. The "more" irrational the value of r, the 
longer the onset of DL is expected to take. In practice 
it is only possible to have rational values of r and thus 
irrational ratios may only be approximated by a close, 
rational ratio [TH ]. This approximation should still give 
similar results to true irrational ratios on the timescale 
investigated in this paper: for a given r = a/b, which is an 
irreducible fraction, the periodicity of the system cannot 
be seen before at least b kicks of the primary pulse-train. 
Thus if the number of kicks from the primary pulse-train 
is not large enough, compared with b, then DL cannot 
be observed and r may be considered to be irrational on 
that timescale. 

We investigate the effect that changing the initial phase 
offset has on DL. We also present evidence of classical 
effects, which in fact dominate the dynamics and are also 
strongly affected by changing the initial phase offset. 

II. ATOM OPTICS KICKED ROTOR WITH 
TWO-FREQUENCY DRIVING 

The Hamiltonian for an AOKR kicked with two driving 
frequencies, corresponding to periods T\ and T2 = rT\, 
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is given, in scaled units, by: 
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where = 2kLX and p = {k/2hkL)p are operators rep- 
resenting the (scaled) atomic position and momentum 
respectively, n\ {n 2 ) is the kicking strength of pulse- 
train 1 (2) (i.e. the area under the pulse), k\ (k 2 ) is 
the pulse shape function for pulse-train 1 (2), ao is the 
scaled time between the first pulse of pulse-train 1 and 
the first pulse of pulse-train 2 (tfj = ao ■ 360° defines 
the initial phase offset) and r = t/T\ is the scaled time. 
We note the scaled commutator relationship [<p, p] = iU, 
where k = Ahk^Ti/m is the effective Planck's constant. 
Note that Xfik^ = k/2 in scaled units. We also define 
M + N = iVtot as the total number of kicks. Note that 
our definition of the ratio is different to that in [j| , where 
r' = hlh = Ti/T 2 = l/r. 

Ntot may be chosen and M and N adjusted such that 
the temporal length of each individual pulse-train is min- 
imised. As pulse overlaps may occur whilst changing the 
phase and the ratio, iVtot may be different from the num- 
ber of individual pulses obtained from the sum of the two 
pulse-trains, N res . 



III. EXPERIMENTAL SETUP 

We realize the AOKR by first trapping and cooling 
approximately fO 5 Caesium atoms to around 5/iK in a 
MOT. The atoms are then subjected to a pulsed stand- 
ing wave of laser light (the kicking beam). The kicking 
beam is obtained by using a f50mW (slave) diode laser 
that is injection locked to a frequency stabilised (master) 
diode laser. The frequency of the master laser may be 
adjusted over a range of about ±2 GHz with respect to 
the D2 line at 852 nm. The output from the slave laser 
passes through an acousto-optic modulator (AOM) that 
is used to produce optical pulses, and this light is coupled 
into an optical fibre and transported to the MOT. The 
kicking laser has a maximum available power of 50mW 
just before the MOT. The standing wave is created by 
retroreflecting the light and has an intensity distribution 
that is well described by a Gaussian with width 0.72mm 
at the position of the cloud. 

An experiment consists of extinguishing the MOT 
beams and quadrupole magnetic field, exposing the 
atoms to a resultant pulse train and allowing the cloud 
to expand for a fixed time (12ms beginning from the mo- 
ment the MOT beams are turned off) so as to increase 
the spatial resolution and to allow for time of flight mea- 
surements. Finally an optical molasses is applied to the 
atoms, "freezing" them in place and causing them to fluo- 
resce, and an image of the spatial distribution of this flu- 



orescence signal is recorded using a CCD camera. From 
this atomic position distribution, and with knowledge of 
the time of flight and the initial position distribution, 
the momentum distribution, and hence the energy of the 
cloud, may be deduced. 

For two- frequency pulse-trains, it is possible for the 
optical pulses to overlap and as such it is important that 
we quantify the exact nature of the pulse train. The 
switching of the AOM is controlled by a programmable 
pulse generator (PPG). The PPG has a timing resolution 
of 40ns and can store a pulse-train of up to 2.6 ms in 
length. Pulse trains of arbitrary shape may be uploaded 
to the PPG by a computer. If an overlap was to occur, 
the maximum laser intensity at the MOT would need to 
be twice as large as for the non-overlapping case. For 
this reason the maximum laser power used was 25mW 
for non-overlapping pulses. 

The pulses output from the PPG were square, with 
negligible rise and fall time. The period of the first 
pulse-train was chosen to be 30//s (k = 3.12) and the 
pulse- width for both trains was set to 480ns. For all ex- 
periments the height of each pulse-train was chosen to be 
identical, and thus n\ = k 2 (although overlaps may cause 
a partial doubling of the height and/or lengthening of a 
resultant pulse). The number of pulses was chosen such 
that Atot = 30. The shape of the pulses at the MOT is 
discussed in section ITVAI 

In order to reduce the uncertainty due to small fluc- 
tuations in the experimental setup (kicking laser power, 
vibrations, slight differences in the initial cloud, etc.), 
three images were taken for each set of parameters. Af- 
ter an experimental run was completed, it was repeated 
and an average was taken over all sets of results, giving a 
measure of the repeatability of the experiment. As a final 
test of repeatability, one set of experiments was entirely 
repeated two weeks later. The results from this second 
set of experiments were identical to the results from the 
first within the experimental error. 



IV. SIMULATIONS 



Both classical and quantum numerical simulations 
were performed for this system, based on the Hamilto- 
nian given by Eq. QJ. The Jacobi elliptic function solu- 
tion for a classical rotator was used to simulate the classi- 
cal system and the quantum system was simulated using 
the Monte Carlo Wavefunction method. In both cases, 
several non-ideal aspects of the physical experiment were 
modelled. These factors were spontaneous emission, the 
spread in kicking strength due to the magnetic sublevels 
of Caesium and the finite sizes of the kicking laser and 
the cloud, and the shape of the optical pulses. 



3 



A. Classical Simulations 

For the classical simulations, free evolution was mod- 
elled by the map: 

<f>[n+l]- = 4>n + + (T[ n+ l]_ ~T n+ )p n+ , (2) 
P[n+1]_ = Pn + , (3) 



r n+ refers to the scaled time between 



where r [n+1] _ 

the end of the n th pulse and the start of the (n + l) th 
pulse. During the time that the potential was on, the 
Jacobi elliptic function solution for the classical rotator 
was evaluated numerically. 

Both the initial position within one potential well and 
the initial momentum of the particle are chosen at ran- 
dom for each trajectory. The position is chosen from a 
uniform distribution in the range [—it, 7r) and the mo- 
mentum is chosen from a Gaussian distribution with a 
standard deviation equal to that given by our measure- 
ment of the temperature. 

The finite size of our cloud, as well as our kicking beam, 
have also been taken into account in the simulation. This 
is done by randomly selecting a position from a Gaussian 
distribution that has the same a as our cloud and evalu- 
ating the intensity at that position. Knowing the relative 
sizes of the beam and the cloud, it is possible to scale the 
kicking strength accordingly. 

Spontaneous emission is taken account of in the sim- 
ulation in the following manner: A number is randomly 
chosen from a uniform distribution between and 1. If 
that number is below our calculated spontaneous emis- 
sion rate per pulse, then a spontaneous emission event 
is deemed to occur during the pulse. The change in mo- 
mentum due to the spontaneous emission is chosen from 
a uniform distribution in the range fifc^) (corre- 

sponding to [— k/2,k~/2) in scaled units). 

It is important that enough timesteps are taken 
throughout a pulse to accurately sample the pulse shape. 
It is also important in the case of overlapping pulses that 
more time steps are taken as the overlaps may cause a 
lengthening of a particular resultant pulse. 

The rise and fall time of the pulses used in experiment 
was also taken account of in simulation. The pulse shape 
was measured just before the MOT, usin g a photodiode 
with a Ins response time. As in reference [l2| , the pulses 
are described by the formula: 
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where fc max is the maximum height of the pulse, t% — t\ 
is the full width at half maximum, St! is the rise time 
(defined such that a straight line going from to 100% 
of the pulse height in time St! matches the slope of the 
rising edge at the half-maximum point), St 2 is the fall 
time (defined similarly to the rise time) and erf (x) is the 
error function. It should be emphasised that, as in |12| . 
this is merely an empirical description of our pulse shape 
and does not arise from physical considerations. 



When the height of the pulse rises to around 10% of its 
maximum, the pulse is deemed to have begun. Similarly 
when the height of the pulse has fallen to less than about 
10%, the pulse is deemed to have ended. Everywhere else 
the pulse height is set to and free evolution is assumed. 

The parameters obtained for our pulses were: <5ti = 
104ns, St2 — 121ns and t%— ti = 396ns. These result from 
the measured power values just before the MOT with the 
PPG outputting a square pulse of duration 480ns. 

As a trade-off between sampling accuracy and CPU 
time, it was decided that at least 16 time steps would 
occur during each non-overlapping pulse. If an over- 
lap occurs so that a pulse is actually longer, then more 
timesteps are chosen such that the length of the timesteps 
are as close as possible to those for the non-overlapping 
case, but are equal or shorter in stepsize. We found, from 
performing simulations, that the energy as a function of 
the number of steps (for a given number of kicks) varied 
only a small amount when the number of time steps was 
10 or greater. 

Implementing this pulse shape, we found that the en- 
ergy output from the simulation was lower compared with 
the result using a square pulse shape. This lower en- 
ergy was in better agreement with the experimentally ob- 
tained energies than the energies obtained using a square 
pulse in the simulation. Also the structures obtained 
(if there were any present to begin with) were slightly 
smoothed out and peaks and dips were broadened (again 
only by a small amount) in comparison with using a 
square pulse of the same area. 

To ensure the area under the real pulse was the same 
as the square pulse, the maximum height of the pulse 
(denoted as fc max in Eq. Q}) was adjusted appropriately. 
For a square pulse the kicking strength is given by k = 
ka, where a is the temporal width of the pulse and k = 
fc(r) is the height of the pulse. When implementing the 
real pulse shape, k still has the same physical significance, 
but must now be expressed in the more general form: 



fc(r)rfr, 



(5) 



which for k(r) = const, reduces to the form given previ- 
ously for the square pulse. 

Due to the fact that we are dealing with a classically 
chaotic process, slight differences in initial conditions give 
wildly different final energies. As such, a very large num- 
ber of trajectories need to be calculated so as to accu- 
rately map out the phase space. In the classical simu- 
lations 10, 000 or 40, 000 trajectories were evolved (more 
trajectories were needed for larger kicking strengths to 
give similar errorbars to the lower kicking strength case). 
This was found to be an acceptable trade-off between 
simulation time and accuracy. 
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B. Quantum Simulations 

For the quantum simulations the method used was sim- 
ilar to that of ref. 13]. A momentum is chosen from a 
Gaussian distribution in the same way as in the classical 
simulation. This momentum is then rounded to the near- 
est momentum ladder state and this is the initial state 
that is evolved. The difference between the momentum of 
the ladder state and that of the chosen initial momentum 
is stored as the quasi-momentum q. 

The initial state is evolved according to the Hamilto- 
nian: 



in between pulses and: 

H = ( ^ g)2 +fc(r)cos(0) 

_ ? £fc(T)^(l + cos(0)), (7) 

where r\ is the probably of a spontaneous emission event 
occuring during one pulse, while the potential is on. Note 
that both k and rj are functions of time (as the intensity 
of the potential changes as a function of time due to the 
pulse shape described by Eq. |0J). The last term in Eq. 
@ exponentially decreases the norm of the wavefunction 
so that spontaneous emission may be included in the sim- 
ulation. When the norm of the wavefunction drops below 
a randomly chosen value, a spontaneous emission event 
occurs during that pulse. 

This process is repeated for multiple initial states (typ- 
ically in the order of 1000) and an incoherent average is 
taken over all of these trajectories to obtain the final mo- 
mentum distribution. 



V. EXPERIMENTAL AND SIMULATION 
RESULTS 

In this section the experimental and simulated results 
will be discussed. These have been divided into three 
sections. First commensurate frequencies in section IV Al 
then incommensurate frequencies in section IV Bl and fi- 
nally the energy as a function of the ratio, r, in section 

fy"Hl 



A. Commensurate Frequencies 

For rational r, giving a periodic pulse-train (which con- 
tains a sufficiently large number of kicks for the period- 
icity to be resolved), coherent quantum effects, such as 
DL, should be observable. 
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FIG. 1: Experimental energies along with quantum and clas- 
sical simulations for r — 1, jc = 17.7, TV-tot = 30 and r\ = 3% 
as a function of ipo. Errorbars are typically ±1.3 energy units 
for the quantum simulation and ±1.5 energy units for the 
classical simulation. 



1. r = 1 

Figure ^ shows experimental results, along with quan- 
tum and classical simulations, for r = 1, Kx = Ki = 17.7 
and Ntot = 30 as the initial phase offset ipo is varied. 

The experimental and simulation results show that 
changing the initial phase offset does have a large effect 
on the resultant energy. A lot of large scale structure, 
as well as various resonance peaks, is seen in the experi- 
mental results and also in the simulations. 

The agreement between the quantum simulation and 
the experimental results is very good. The energies from 
the classical simulations are larger than those of the 
quantum simulations for most values of ipo , suggesting 
the presence of DL. 

Most of the large scale structure may be explained by 
the presence of KAM boundaries and the overlapping of 
pulses. As -00 tends to moderately small or moderately 
large values, the resultant pulse-train consists essentially 
of double pulses (where the spacing between each pair 
of pulses is much larger than the spacing between each 
pulse in the pair). It is well known that double pulses 
produce KAM boundaries (which are a classical effect), 
and this situation has previously been used to study the 
effects of quantum tunnelling through KAM boundaries 
[TH Il6| . KAM boundaries should therefore explain the 
large scale decrease in energy in these regions, even for 
the quantum simulation. This argument is supported by 
viewing the momentum distributions in these regions (as 
an example see figure 0). KAM boundary effects seem to 
become less important in the region 75° ~< ipo ~< 285° 
where the energy seems to take on a more uniform value 
(ignoring resonance peaks). 

The large peaks as ijjQ - ► 0°, 360° are explained by the 
overlap of pulses. At ipo = 0°,360° the pulses are over- 
lapping completely (periodic with period 30/^s) and there 
are only half as many resultant pulses as there would have 
been without overlaps. Also k is twice as large as with- 
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FIG. 2: Experimental KAM boundary momentum distribu- 
tion for r = 1, Vo = 34.3°, k = 10.1 and AT tot = 30. The 
shoulders are a symptom of KAM boundaries and appear at 
the momentum values that KAM boundaries occupy. Trajec- 
tories initially inside the boundary (essentially all trajectories 
are inside for a cloud temperature of around 5/iK) have dif- 
ficulty diffusing through it, hence the fairly sharp drop off in 
population. Note that this distribution has not been decon- 
volved with the initial position distribution, but this has only 
a slight effect on its shape. 



out overlaps. Ignoring anomalous diffusion (and DL), the 
diffusion rate increases as n 2 and the energy increases 
roughly linearly with kick number. So with twice the 
kicking strength and half the number of pulses, the final 
energy obtained would be roughly two times as large as 
if there were no overlapping pulses. This rough calcu- 
lation is in approximate agreement with what is seen in 
our results. 

It is interesting to note the positions of the resonance 
peaks (e.g. at ip = 90°, 120°, 180° in figure [TJ. They 
occur at rational fractions of 360°, i.e. when ao (from 
Eq. is a rational fraction. These observations are 

suggestive of mode locking 0] . The larger peaks appear 
to follow the general pattern of the Farey tree [Tlj in 
that the simpler the rational fraction, the larger the peak 
seems to be. 

Figure shows the results with k\ — K2 — 10.1 and 
all other parameters the same as previously. A similar 
large scale structure is obtained to that when k = 17.7 
- large peaks at ipQ — * 0° , 360° , KAM boundary regions 
for moderately small and large ipo and resonance peaks. 
There are some differences however. Some of the reso- 
nance peaks from figure ^ are now anti-resonance dips 
(e.g. tpo — 90°, 180°) and there is a broad peak (with a 
dip in the middle) centred on ipo = 180°. 

The distinction between quantum and classical be- 
haviour is not so clear in the experimental data as for 
k = 17.7. However, similarity with the quantum sim- 
ulation is observed in the regions ipa — > 0°,360° and 
120° < tpo < 240°. The qualitative differences between 
the quantum and classical simulations are expected to be 
largest in these regions, as the pulse train becomes more 




FIG. 3: Experimental and simulated (both quantum and clas- 
sical) energies for r = 1, k = 10.1, Atot = 30 and rj = 2.8% 
as a function of ipo- Errorbars are typically ±1.2 energy units 
for the quantum simulation and ±1.7 energy units for the 
classical simulation. 



periodic for these initial phase offsets. 
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FIG. 4: Classical simulation for k = 10.1 and -q = 2.8% as a 
function of ipo for r — 1 and 100 kicks. Errorbars are typically 
2.3 energy units for the classical simulation and ±3.7 energy 
units for the quantum simulation. 

Quantum and classical simulations were run for k = 
10.1 and iVtot = 100 (figure to see how these 
peaks/dips evolve in the quantum and classical cases. In 
the classical simulation all peaks/dips are present and 
become sharper. In the quantum simulation, however, 
only the peaks/dips that were largest for A^ot = 30 are 
present - the other resonances seem to have been washed 
out. This suggests that the origin of these peaks/dips 
is classical, not quantum mechanical, and that quan- 
tum effects may destroy these resonances. Another point 
to note is that the relative difference in energy between 
the quantum and classical simulations is much larger for 
Ntot — 100 than it was for Ar to t = 30, suggesting the 
presence of effects of DL. Note that atomic spontaneous 
emission leads to a finite late-time diffusion rate in the 
quantum case, but this rate is generally slower than the 
corresponding classical rate. Note also that, despite this 
breakdown of DL, the momentum distributions can re- 
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main essentially exponential in shape even after this num- 
ber of kicks [13 • 

In the classical simulation the peaks/dips are clearer 
after more kicks. This is to be expected as the periodicity 
is resolved better after more kicks. Furthermore, figure 
0] shows that if ceo is written as an irreducible fraction, 
czq = a/b, the peaks always appear when b is an odd 
number and the dips always appear when b is an even 
number - regardless of whether a is even or odd. This 
observation held true even for a classical simulation for 
Ntot = 500, where peaks corresponding to less rational 
values of ceo were able to be resolved above the noise. 

Poincare sections were plotted to search for structures 
that may explain the resonances/ anti- resonances. Al- 
though there were structures present, clear correlation 
with the shape of the energy curve was not obvious. Gen- 
erally the island structures in phase space remained rea- 
sonably constant over large ranges of tpo values. An ex- 
planation for both the appearance and positions of these 
peaks and dips will require further investigation. 

2. Momentum Distributions for r = 1 

As mentioned earlier, the difference in energy of the 
classical and quantum simulations suggest the presence 
of DL for both values of k. The key feature that distin- 
guishes DL from other effects, however, is an exponen- 
tially localised momentum distribution. Figure [5] shows 
the momentum distributions obtained for k — 10.1 and 
tpo = 180°. It still shows an essentially exponential mo- 
mentum line shape, even for the level of spontaneous 
emission present, in good agreement with the quantum 
simulation. A Gaussian distribution is seen for the clas- 
sical simulation. The distributions showing DL are most 
prominent around tp = 120°, 180°, 240°, but are seen to 
a lesser degree for almost all values of tpo- 
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FIG. 5: DL momentum distribution for r = 1, k = 10.1, 
Ntot = 30, T) — 2.8% and tpo — 180°. Note the experimental 
distribution has not been deconvolved with the initial position 
distribution, but this has only a slight effect on its shape. 

In the region 110° < ipo < 250° there is a mix- 



ture of DL and KAM boundary distributions (see fig- 
urc[S]where an exponential distribution with shoulders is 
found, which is significantly narrower than the classical 
distribution and in agreement with the quantum distri- 
bution), except at tpo = 120°, 180°, 240° where only DL 
is seen. This suggests that KAM boundaries may even be 
having an effect on the dynamics in this region (it should 
be noted that the diffusion rate through a cantorus can 
be lower quantum mechanically than classically p^lll6j|). 
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FIG. 6: Mixed momentum distribution for r = 1, k = 10.1, 
Ntot = 30, r\ — 2.8% and ipo = 170°. Note the experimental 
distribution has not been deconvolved with the initial position 
distribution, but this has only a slight effect on its shape. 

Interestingly, at ip = 120°, 180°, 240° the resultant 
pulse-train contains (or very nearly contains) completely 
periodic pulses (i.e. having equal spacing). The re- 
sultant pulse-train for tpo = 120° may be viewed as 
being completely periodic with period T — 10/is, but 
with every third pulse missing. The resultant pulse-train 
for tpo — 180° is completely periodic with T = 15/xs. 
tpo = 240° is a similar case to tpo = 120°. 

The largest qualitative difference between the quantum 
and classical simulations is around tpo = 180°. Although 
in figure|3the qualitative structure is very similar for the 
two simulations, in figure nt nere is a dip at tp = 180° 
for the quantum case and only a peak in the classical 
simulation. This region was looked at in more detail for 
various kicking strengths. The results are shown in fig- 
urc0and indicate a reasonably complicated structure for 
different values of k for both quantum and classical dy- 
namics. The energy is lower in the quantum simulations 
(the peak has a dip in it), compared with the classical 
simulations, as tpo — > 180°. This can only be explained 
by DL and indeed the momentum distribution exhibited 
an exponential lineshape. Observing strong DL at this 
value of the initial phase is not surprising because, as 
mentioned earlier, the pulse-train is completely periodic. 
Experimental results are always in qualitative agreement 
with the quantum simulation and show DL particularly 
around tp = 180°. 
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FIG. 7: Results for r = 1 and TVtot = 30. Classical sim- 
ulations have the thickest lines, quantum simulations have 
thinner lines and experimental results have the thinnest lines 
with errorbars. 



3. r = 1.4 
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FIG. 8: Experimental energies along with quantum and classi- 
cal simulations for r = 1.4, n — 10.1, Mot = 30 and rj = 2.8% 
as a function of ipo. Errorbars are typically ±1.7 for the quan- 
tum simulation and ±1.9 for the classical simulation. 



Results for other two-frequency ratios were also ob- 
tained. The experimental results for r = 1.4 with 
k = 10.1 and k = 17.7, along with simulation results, are 
shown in figures [S] and |§1 respectively. Again the quali- 
tative agreement between experiment and simulations is 
good. There is definite structure present for both kicking 




FIG. 9: Experimental energies along with quantum and clas- 
sical simulations for r — 1.4, k — 17.7, iVtot = 30 and rj — 3% 
as a function of ipo for. Errorbars are typically ±1.7 energy 
units for the quantum simulation and ±2.6 energy units for 
the classical simulation. The experimental points have been 
joined to show the similarity in structure with the simulated 
results. 



strenghs which is shown in the experimental result as well 
as simulations. However, for this value of the ratio there 
are no noticeable resonance peaks or anti resonance dips. 
Also, there is no difference in energy between the classi- 
cal and quantum simulations within numerical error, ft 
seems that for r = 1.4, 30 kicks is not enough to resolve 
the periodicity sufficiently for DL to have an effect. This 
is backed up by the momentum distributions which do 
not show any evidence of DL. They do, however, show 
that KAM boundaries are present. 

These results show that even for this reasonably com- 
mensurate frequency, the periodicity is not resolved well 
enough on this time scale for the system to exhibit sig- 
nificant quantum effects. 



B. Incommensurate Frequencies 

If the ratio of pulsing periods is an irrational number, 
DL is not expected as the accumulated quantum phase 
from kick to kick will be approximately incoherent (for 
a sufficiently large number of kicks). Thus the dynamics 
should be dominated by classical behaviour. 

Figures 1101 and 1111 show simulated and experimental 
results for r = \[2 with n = 10.1 and k — 17.7, respec- 
tively, after iV to t = 30 kicks. Both quantum and classical 
simulations are shown. 

As expected the energy shows a more uniform value 
compared to rational ratios. Contrary to what may have 
been expected, the energy does show some dependence on 
rpo and is not simply the flat line expected in the zero cor- 
relation limit. There is definite structure present for both 
classical and quantum simulations, and the experimental 
results show this same structure. The classical simula- 
tion is more noisy but despite this, the broader features 
of the quantum and classical simulations are identical 
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FIG. 10: Experimental energies along with quantum and clas- 
sical simulations for r = \/2, k — 10.1, jV to t = 30 and 
77 = 2.8% as a function of ipo. Errorbars are typically ±1.7 en- 
ergy units for the quantum simulation and ±2.1 energy units 
for the classical simulation. 
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FIG. 11: Experimental energies along with quantum and clas- 
sical simulations for r — \/2, K = 17.7, N to t = 30 and r\ — 3% 
as a function of tfio- Errorbars are typically ±1.6 energy units 
for the quantum simulation and ±2.6 energy units for the clas- 
sical simulation. The experimental points have been joined to 
show the similarity in structure with the simulated results. 



within error. This similarity could not be clearer than in 
the higher kicking strength case (see inset in figure HT|) . 

The variations in energy, for both cases, seem to have 
some sort of systematic variation with initial phase. 
When an FFT is performed the main spectral compo- 
nent is found to be l/A-00 ~ 1/70° for k = 10.1 and 
l/AV>o ~ 1/30° for k = 17.7. The variations in energy, 
while small compared to the total energy (in the order 
of 10%), are nonetheless present. The mean energy is 
around 41 two photon recoils for n = 10.1. This corre- 
sponds to the zero-correlation (quasilinear) limit, taking 
into account the spread in kicking strengths present in 
the real experiment due to the finite size of the beam 
and the cloud, as well as the magnetic sublevels of Cae- 
sium, and also the pulse-shape. 

Figure^2 a l so shows that the energy of the experimen- 
tal data is lower than both the quantum and classical 
simulations. This suggests that, at least for higher kick- 



ing strengths, we may be losing some of the atoms in the 
kicked cloud during the molasses freezing phase before 
imaging. This energy offset can also be seen in other 
figures with higher kicking strengths. 

The energy variations may be due to transient early 
time effects that have not yet been "washed out" after 
30 kicks. Only a finite length pulse train can be applied 
to the rotors, and in theory an infinitely long train would 
be needed to realise an incommensurate frequency ratio 
properly. 

This is related to the way in which an irrational num- 
ber may be approximated arbitrarily well by an infinite 
sum of rational fractions. The number of terms in the 
sum determines how well the irrational number is ap- 
proximated. The continued fraction representation of 
any number, u>, between and 1 may be written as: 



An irrational number will have infinitely many non-zero 
integer terms, en. The more terms that are included in 
the continued fraction, the better the approximation to 
the irrational number. If an integer offset, ao, is included, 
then any irrational number may be approximated by this 
method. 

In terms of the aforementioned pulse train, due to the 
discretisation of time, more pulses is analogous to having 
more ai terms in the continued fraction representation. 

Figure IT21 show results for r = \fi with all other pa- 
rameters the same. Slightly different structure is seen to 
that for r = \/2, however the same trends are exhibited as 
for r — y/2, i.e. periodic variations on an approximately 
uniform background energy. 
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FIG. 12: Experimental energies along with quantum and clas- 
sical simulations for r = \/3, k = 10.1, N to t = 30 and 
r\ — 2.8% as a function of ipo. Errorbars are typically ±1.6 en- 
ergy units for the quantum simulation and ±2.1 energy units 
for the classical simulation. 

If the cause of the structure seen is the 'small' number 
of kicks in the pulse-train, then having more kicks should 
cause the variations to be washed out. It is not possible 
with the current experimental setup to test this theory, 
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as not enough pulses can be stored in the memory of the 
PPG. However it is possible to look at the effects found 
in simulation by increasing the kick number. Figure ITTA 
shows the classical and quantum simulations for k — 17.7, 
r = \/2 and 90 kicks. The structure is clearly different 
from that after 30 kicks (figure EH- Four times the num- 
ber of trajectories have been used here so as to keep the 
errors similar to the lower kicking number. The abso- 
lute height of the variations is approximately the same 
as for the lower kicking number, however the variation 
relative to the total energy has decreased to around 3% 
(compared to 10% for the 30 kick case). There is less 
structure than in the low kick number case and it seems 
as though the energy is tending more towards a more uni- 
form value as the kick number is increased. As more kicks 
are added to the pulse-train, it is likely these transient 
effects would become less and less noticeable. 
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FIG. 14: Simulations of the system described in with 
iVtot = 50 and n — 10. The number of zero velocity atoms 
has been calculated by integrating the normalised final mo- 
mentum distribution over — e to e for some fixed value of e. 



FIG. 13: Quantum and classical simulations of the energy 
for TVtot = 90, r — y/2 and k = 17.7. Errorbars are typi- 
cally ±5.2 energy units for the quantum simulation and ±6.0 
energy units for the classical simulation. 

Weak KAM boundary distributions are seen for prac- 
tically all irrational values of r and all ?po i n the simu- 
lations (however the boundaries are not nearly as strong 
as in the cases with rational r) and no DL is seen in the 
quantum simulation or the experimental results. How- 
ever, the KAM boundaries were too weak to be observed 
in the experimental momentum distributions. It is not 
clear that the KAM boundaries are responsible for the 
observed structure, however it is clear that the dynamics 
for irrational ratios is again governed by classical effects. 

C. Energy Growth as a Function of the Ratio of 
the Driving Frequencies 

As mentioned earlier, when the ratio of the driving 
frequencies was highly rational, Ringot et al. found 
an increase in the number of zero velocity atoms com- 
pared to when it was highly irrational. This behaviour 
was attributed solely to DL. Although the parameters 
of ratio and initial phase have been defined differently 
here than in 9] , they are simply related by r' = 1/r and 



ip'o = V^o/ 7 "- Using these conversions, it can be seen from 
figure El that for r = r' = E, k « 10 and Vo = V'o = 52° 
(the specific phase used in KAM boundaries play a 
large role in the dynamics and in fact also halt the energy 
growth. Hence the number of zero velocity atoms would 
be larger than if KAM boundaries were not present. The 
case where KAM boundaries are weak can be seen in fig- 
ureQKr' = 1/V2 « 0.707) with ip = 52° x V2 « 73.5°. 
Here the energy is larger than the previous case and hence 
fewer zero velocity atoms are expected. This is true for 
both the quantum and classical cases and hence both 
cases should show an increased number of zero velocity 
atoms for r' = 1 compared to r' — y/2. 

To verify this we also ran quantum and classical sim- 
ulations for parameters equivalent to those stated in 0|. 
The results are shown in figure ITlk . In our simulations 
both the quantum and classical cases display peaks in the 
number of zero velocity atoms for rational values of r'. 
Also, we examined the momentum distributions to deter- 
mine whether DL or KAM boundaries were responsible 
for these peaks. 

In contrast to the explanation given in jj], by exami- 
nation of the simulated momentum distributions, we find 
that KAM boundaries, and not DL, are primarily respon- 
sible for the increase in the number of zero velocity atoms. 
Only a small amount of DL is observed for the most ra- 
tional ratios and even then the most dominant effect is 
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that of KAM boundaries. 

The fact that the effects of DL seem to be easily de- 
stroyed by changing the initial phase offset suggests that 
the results obtained in were to some degree due to the 
particular choice of initial phase and may not be true in 
general. To check for this we also ran simulations with 
tp' = 150°. The results are shown in figure l*HTt>. There 
is a difference between this case and the ip' = 52° case 
in terms of structure - particularly around r' — 2. Peaks 
still occur at rational values of the ratio in both the clas- 
sical and quantum cases. The peaks in the classical and 
quantum cases are approximately the same height, ex- 
cept for the peak at r' = 1, which is noticeably lower in 
the classical case. The quantum momentum distribution 
at this value of r' distinctly shows DL, as well as KAM 
boundaries. 



VI. CONCLUSION 

To summarise, for rational ratios, those values of the 
initial phase that give a more periodic resultant pulse- 



train do show DL on the timescale investigated. How- 
ever, the values of the initial phase that resulted in less 
periodic pulse-trains do not show DL on the timescale 
investigated, even for the most rational ratios. 

The initial phase offset of the resulting pulse-train has 
a large effect on the dynamics in both the quantum and 
classical cases for rational ratios, and this effect is en- 
hanced on longer timescales. The only effect present for 
irrational ratios is a weaker transient structure that has 
a diminished effect on longer timescales. 

No DL is seen when the ratio, r, is irrational for any 
value of the initial phase offset. That is, quasi-periodic 
driving destroys quantum kick to kick correlations on the 
timescale investigated. 

Most of the observed effects have been found to be 
classical, not quantum mechanical, in origin and are not 
due to dynamical localisation. A dominant factor in the 
dynamics is KAM boundaries. 
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